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Attraction/repulsion switching of non-equilibrium depletion interaction 
caused by blockade effect in gas of interacting particles. II 
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The effect of concentration-dependent switching of the non-equilibrium depletion interaction be¬ 
tween obstacles in a gas flow of interacting Brownian particles is presented. When increasing bath 
fraction exceeds half-filling, the wake-mediated interaction between obstacles switches from effective 
attraction to repulsion or vice-versa, depending on the mutual alignment of obstacles with respect 
to the gas flow. It is shown that for an ensemble of small and widely separated obstacles the dis¬ 
sipative interaction takes the form of induced dipole-dipole interaction governed by an anisotropic 
screened Coulomb-like potential. This allows one to give a qualitative picture of the interaction 
between obstacles and explain switching effect as a result of changes of anisotropy direction. The 
non-linear blockade effect is shown to be essential near closely located obstacles, that manifests itself 
in additional screening of the gas flow and generation of a pronounced step-like profile of gas density 
distribution. It is established that behavior of the magnitude of dissipative effective interaction is, 
generally, non-monotonic in relation to both the bath fraction and the external driving field. It 
has characteristic peaks corresponding to the situation when the common density “coat” formed 
around the obstacles is most pronounced. The possibility of the dissipative pairing effect is briefly 
discussed. All the results are obtained within the classical lattice-gas model. 

PACS numbers: 05.40.-a, 51.10.+y, 68.43.Jk 


I. INTRODUCTION 

Motion of inclusions or probe-particles through a 
medium is accompanied by the medium perturbation 
(e.g., perturbation of its density) that can manifest it¬ 
self in the form of wakes. The medium perturbation 
can, in turn, induce a non-equilibrium interaction be¬ 
tween the inclusions. Such interaction is responsible, in 
particular, for the coherent part of the collective fric¬ 
tion force as well as for possible formation of dissipative 
structures in an ensemble of the inclusions. The nature 
of the medium perturbation and the properties of the 
induced non-equilibrium interaction are defined by the 
properties of the medium (e.g., by its nonlinearity) and 
the mechanism of energy losses. The perturbation can 
lead to generation of vortices, Cherenkov radiation, or 
local phase transitions (some more effect can be found in 
hydrodynamics M, optics (?!, plasma physics @-0, 
quantum liquids and Bose condensates [l4j-[2C|). 

In dissipative media, the induced non-equilibrium in¬ 
teraction between inclusions can conditionally be divided 
into reactive and dissipative parts. In the simplest case, 
when the speed of a probe-particle is rather small (e.g., 
smaller than the speed of sound in a medium and the 
hydrodynamic effects can be neglected, the medium per¬ 
turbation can be described in the diffusive approximation 
[2ll |. The diffusive wake may be of large spatial and tem¬ 
poral extensions with power-law damping (see [l^-[26[), 
which is an evidence of long-time memory of the medium 
about the particle passage. The long-living wakes of in¬ 
dividual particles lead, in turn, to a long-range effective 
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dissipative interaction between the particles [27j. This 
can be qualitatively described using the linear response 
approximation [lfl . |28I | . However, the linear response ap¬ 
proximation, giving a qualitative picture of medium per¬ 
turbation, leads to incorrect results for wakes, dissipa¬ 
tive interaction and, in general, does not give adequate 
description of non-linear media [20l . 


In the present paper, we will be interested in the dis¬ 
sipative interaction between inclusions induced by their 
wakes in a nonlinear medium, resorting to an example 
of a Brownian gas with short-range inter-particle repul¬ 
sion (the hard-core interaction). In this case, the dissi¬ 
pative interaction between inclusions is often called the 
non-equilibrium depletion or entropic interaction (e.g., 

0 S3, S3). At equilibrium, the depletion interaction 

is usually short-range; its spatial range is of the order 
of the characteristic length scale of the medium parti¬ 
cles [13, 0. In contrast, the non-equilibrium forces be¬ 
tween impurities may exhibit long-range behavior due to 
a long-living diffusive wake induced by their motion ms 
In addition, such forces often have un¬ 


usual properties, e.g., they violate the Newton’s third law 
S3, S3, S3, Ss] . The non-Newtonian behavior of the non¬ 
equilibrium depletion force was demonstrated at low gas 
concentrations [27!], when interaction between gas parti¬ 
cles is negligible. 


To describe the non-equilibrium depletion force for a 
gas of interacting particles, we turn to the simplest model 
of a lattice gas, when each lattice site can be occupied 
by only one particle. Even such a short-range interac¬ 
tion results in a number of unexpected kinetic effects, 
e.g., the “back correlations” effect j56|, drifting spatial 
structures S3 - l39j] . effects of “negative” mass transport 
[13-0, induced long-time correlations [43], and the dis¬ 
sipative pairing effect for tracers passing through a lattice 
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gas [44|. Increasing of gas concentration (bath fraction) 
leads to enhancement of the role_of interaction between 
gas particles. As was shown in (45!), this implicates sig¬ 
nificant changes in the shape of wake of a fixed inclusion 
(or an obstacle) in a gas flow — wake inversion. In turn, 
the wake-mediated interaction between obstacles should 
be sensitive to the crucial transformation of the wake 
structure. 

In this paper we examine how the short-range repul¬ 
sive interaction of gas particles affects the behavior of 
dissipative forces between obstacles embedded into the 
gas flow. In particular, we show that increasing of gas 
concentration can lead to the switching, or sign change, 
of the effective dissipative interaction between obstacles 
to its opposite, e.g., from attraction to repulsion or visa- 
versa. This effect is entailed by the obstacle wake in¬ 
version considered in [45j. For closely located obstacles, 
the interaction of gas particles is shown to provoke non¬ 
linear blockade effect, that results in formation of a com¬ 
mon “coat” of gas density perturbation around the ob¬ 
stacles with a pronounced step-like behavior of its distri¬ 
bution. In turn, the common non-linear coat can signify 
the dissipative pairing of the obstacles, see [44|. In the 
case of small and widely separated inclusions, when the 
non-linear effects are less significant, we show that the 
dissipative interaction between them belongs to the type 
of induced dipole-dipole (generally multipole) interaction 
associated with anisotropic screened Coulomb-like poten¬ 
tial. To demonstrate the above mentioned phenomena, 
we use the mean-field and the long-wavelength approxi¬ 
mations, neglecting the short-range correlations and fluc¬ 
tuations in the gas, see [ill |45|. 

Our paper is organized as follows: In Sec.|TI]we specify 
the kinetic equations to be used and briefly discuss the 
employed approximations. The main results on dissipa¬ 
tive forces are contained in Sec. m In Sec. 1III A1 the case 
of small (point-like) inclusions is considered in the linear 
flow approximation. In Sec. IIII51 the non-linear blockade 
effect (i.e., screening of gas flow) is discussed for large and 
closely located obstacles. In Sec. IIII Cl the case of two 
moderately separated obstacles is considered numerically 
for two spatial configurations. Sec. IIVI briefly summa¬ 
rizes obtained results. Appendixes contain the outlines 
of two analytic approaches used in Sec. IIII A] a naive one 
(Appendix [A} , giving a rough sketch of the dissipative 
interaction behavior, and a more sophisticated one (Ap- 
pendix[B]), based on the single-layer potential method for 
inclusions with sharp boundaries. 


II. MODEL 

As was shown in [HI j4f|, an obstacle in a lattice 
gas flow can be considered as a limiting case of a two- 
component gas: one of the components is static while 
the other one is mobile and driven by a uniform external 
field. We employ the simplest model of a two-component 
lattice gas, when each lattice site can be occupied by 


only one particle, see f§6| . Kinetics of a multicomponent 
lattice gas is defined by the jumps of its particles to the 
neighboring vacant sites. The variation of the ith site 
occupancy by the particles of the ath sort during the 
time interval At, To <C At tj (to is the duration of a 
particle jump to a neighboring site and r; being the life¬ 
time of a particle on a site), is described by the standard 
continuity equation (see, e.g., [§§1 |46}~) 

<(t + At) - <(t) = £ (J“ - J“) + 6J?, (1) 

j 

where a and /? label the particle species and nf = 0,1 are 
the local occupation numbers of the ath particles at the 
tth site. J“- = zAnf ^1 — At gives the average 

number of jumps of the ath particles from site * to a 
neighboring site j per time At. zA = v a is the mean 
frequency of these jumps. The term <5J“ = ~ 

(5J°j) stands for the Langevin source that is defined by the 
fluctuations SJ°i of the number of jumps between sites j 
and i during At j46|. These fluctuations are caused by 
fast, compared to the time scale At, processes and will 
be neglected for simplicity. It means that we disregard 
the fluctuation-induced forces. 

In what follows we consider only two components, mo¬ 
bile and static, which are labeled by n and «, respectively. 
In the absence of external fields we suggest for a regular 
lattice that zA = v = const for the component n, while 
the component u is assumed to be at rest, /A = 0. The 
presence of a driving field leads to asymmetry of the par¬ 
ticle jumps. Assuming the activation mechanism of the 
jumps and a weak driving field G, frequency Uji may be 
written as zA ft ^[1 + (G, r* — rj)/(2fcT)], or v± ft: v±8v, 
where u + and v~ denote the jump frequencies along and 
against the field, respectively. 8v = vf\G\/ (2 kT) (i is the 
lattice constant), condition tjG|/(2 kT) < 1 is assumed 
to be satisfied. 

Equations for the average local occupation numbers 
can be obtained from Eqs. (JT|) using the local equilibrium 
approximation (the Zubarev approach) [46j, [47} which 
coincides, in our case, with the mean-field approxima¬ 
tion [48} . Introducing time derivatives [4^ , in the long- 
wavelength approximation (see f§i] - l39[ IT 11 1 the macro¬ 
scopic kinetics of the mobile component n is given by the 
equation 

d T n = V 2 n —V(izVn —nVzz) — (g, V)[?r(l — u — n)], (2) 

where n = n( r, r) and u = u( r) are the average occu¬ 
pation numbers of the two components at the point r 
(0 < n < 1 and 0 < u < 1) and g = tG/(2kT). Here, 
we have introduced the dimensionless spatial coordinate 
r/(. —> r and time r = vt, and d T stands for the partial 
time derivative. Note that equations of the form ©, as 
well as their generalizations for two- and multicomponent 
systems, also appear in the problems of nonlinear cross¬ 
diffusion with size exclusion [50|, diffusion in monolayers 
of reagents on the surface of a catalyst [5l| and serve as 
a model of fast ionic conductors [§71] ■ 
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In the non-equilibrium case, there are various ap¬ 
proaches to introduce the dissipative force (or interac¬ 
tion) between inclusions via the Brownian gas environ¬ 
ment. The approaches are not equivalent to each other 
and may lead to different results in general, see [3fl]. To 
introduce the force acting on an obstacle, we first con¬ 
sider a point-like inclusion (impurity) occupying a lattice 
site Rj with a given interaction potential t/(r, — Rj) be¬ 
tween the inclusion and a particle of the lattice gas at site 
Then, Hamiltonian of the lattice gas in the presence 
of impurities is written as H = H 0 + H int , where H 0 is 
the Hamiltonian of the lattice gas without inclusions, and 
Hint = — Rj) describes interaction between 

gas particles and impurities, m = 0,1 is the occupation 
number of site rj. 

At equilibrium, the total force acting on the j th inclu¬ 
sion can be written as (see 0) 

= {~SRj Hint ) = 2 P({n}) (3) 

{n} 



where p({n}) is the equilibrium probability (or statisti¬ 
cal operator in the matrix representation 1461 ]) of a given 
occupancy configuration {n} 

p({n}, 0) = Z^ 1 exp(— H{n}/kT) (5) 

and Z = ]P{ n } exp(— H{n}/kT); (rij) is the mean occu¬ 
pation number at site rj that describes the equilibrium 
distribution of gas concentration. The force /J 9 , Eqs. © 
or 0. can be expressed in terms of the gas free energy 
F = — kT In Z as 

F = -m~ F ' <6) 

This relation is often used to define the equilibrium de¬ 
pletion force jH, [53]. 

In this paper, we use another approach based on ex¬ 
pression © written with non-equilibrium statistical op¬ 
erator pt({n}) (see HI) 

F” = E <’<((’*» = Y.^>£- U ( r - - r j>. 

{n} x J ' i 

( 7 ) 

where pt({n}\ obeys a master equation for the hopping 
process, see HI], and (rij) t = J^p n yriip t ({n}) is non¬ 
equilibrium gas concentration. Yet another approach 
consists in generalizing Eq. © to the non-equilibrium 
case by introducing an effective non-equilibrium poten¬ 
tial or non-equilibrium free energy for a gas (30;, [13, IH3 - 
[56} . As was shown in [30], these two definitions of the 
non-equilibrium force are not equivalent. Representation 
© for the force exerted by gas particles on an obsta¬ 
cle is similar to the hydrodynamic definition of the force 
which, in particular, was used in [27j to describe the non¬ 
equilibrium depletion interaction between obstacles in a 


gas of non-interacting particles. Here, we use representa¬ 
tion © to describe the non-equilibrium depletion forces 
acting between obstacles via gas perturbation. 

In the continuum limit and the mean-field approxima¬ 
tion, /" e? takes the form 

/" e9 = - J U(r — Rj) V r n(r, t) dr, (8) 

where n(r,t) = (n(r)) t . When the obstacle is a cluster 
formed by particles of the second (heavy) gas component, 
potential U(r) describes the concentration distribution 
of that component and n(r,t) obeys Eq. © obtained 
in the long-wavelength approximation. In what follows, 
to separate out the contribution of the gas perturbation 
5n(r, t) induced by the gas flow (or the external field g) 
from the total force ©, we consider the quantity 

fj=-Ju(r-Rj)V r 6n(r,t)dr, (9) 

where Sn(r,t) = n(r,t) — no(r), no(r) is the equilibrium 
concentration distribution, and no(r —> oo) —> no = 
const stands for the average equilibrium concentration of 
gas (fraction of the full lattice occupation, 0 < no < 1). 

In the case of inclusion with a sharp boundary, fj takes 
the conventional form 

fj = ~ [ n(r)5n(r) dr, (10) 

Js d 

where Sj is the surface of jth inclusion and n(r) is its 
exterior normal at the point r. In what follows, we will 
be interested in non-equilibrium steady-state interaction, 
i.e., in the limiting case t —> oo. We will use the lattice 
gas model © in the mean-field approximation (neglect¬ 
ing the fluctuation part) and its continuum version © 
to describe the character of the dissipative interaction 
between obstacles. 


III. INVERSION OF WAKE AND SWITCHING 
OF DISSIPATIVE INTERACTION 

In this section we consider how the gas particle in¬ 
teraction and non-linear screening of gas flow affect the 
behavior of dissipative forces acting on obstacles. In par¬ 
ticular, we show that the short-range repulsive interac¬ 
tion between gas particles can lead to switching of non¬ 
equilibrium depletion interaction between obstacles, e.g., 
from effective repulsion to attraction, as the equilibrium 
gas concentration n o increases. Such interaction switch¬ 
ing is directly related to the obstacle wake inversion ef¬ 
fect considered in [45j]. As was shown in (4f|, increasing 
of gas concentration no leads to drastic transformation of 
the inclusion wake structure: typical wake Binma 
with an extended depleted region behind the inclusion 
and localized dense region in front of it [at no < 1/2, see 
Fig. ©a)], acquires an unusual inverted structure with 
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Figure 1. Transverse alignment. Steady-state concentration distributions (average occupation numbers) of the gas particles 
n(x, y) near the obstacles, evaluated numerically within the mean-field approximation of Eq. ([TJ in the 2D case, correspond to 
various equilibrium concentrations no = 0.2 (a), 0.5 (b), and 0.8 (c). Panels (a) to (c) illustrate three different regimes of the 
dissipative interaction: (a) effective repulsion (/f 2 > 0, / 21 < 0, |/ 21 | = |// 2 |), (b) no interaction (/ 21 = /f 2 = 0), (c) effective 
attraction (/ 21 > // 2 ). The external held g (|g| = 0.5) is directed along the z-axis; the impermeable (u = 1) circular obstacles 
are of radius a = 7 (in units of £); the distance between the obstacle centers equals 4a. The gray background corresponds to 
the equilibrium gas concentration no for every contour plot, in consistence with the color bars. Spatial coordinates are in units 
of l. 


an extended dense region ahead of the obstacle and a lo¬ 
calized depleted region behind [at no > 1/2, Fig. [2(c)]. 
Note that the possibility of wake and force switching can 
be easily shown by using the “hole-particle” symmetry of 

Eq. ©, see [410- 

Switching of the wake “direction” at high gas concen¬ 
tration is caused by the enhancement of the role of in¬ 
teraction between gas particles, in particular, can lead 
to the non-linear blockade effect. This effect is signifi¬ 
cant near the obstacle surface, especially for large and 
for closely located ones. For a relatively large obstacle 
and sufficiently high concentration no, the gas flow gen¬ 
erates a dense region ahead of the obstacle as the gas 
particles have no time to leave this zone via lateral diffu¬ 
sion. Such a strong accumulation of the gas particles lo¬ 
cally enhances the significance of the interaction between 
them, so that the dense region ahead of the obstacle has 
to grow. Similar behavior arises for closely located obsta¬ 
cles when their individual density perturbation “coats” 
overlap leading to formation of a common “coat” around 
them and to additional screening of the gas flow. The 
latter means that peculiarities of dissipative interaction 
between closely located obstacles are determined by the 
non-linear blockade effect for which the term ~ n 2 in 
Eq. @ is responsible. We consider these non-linear ef¬ 
fects numerically on the basis of mean-field version of 
Eq. ©, neglecting gas fluctuations. 

In the particular case of relatively small and distant ob¬ 
stacles, the interaction between gas particles can be taken 
into account in the linear approximation [4f|. That ap¬ 
proximation allows one to obtain analytical expressions 
for the asymptotic behavior of both the density pertur¬ 
bation far from obstacles and the dissipative interaction 
between them. Now we proceed to this case in the sub¬ 
section below. 


A. The asymptotic behavior of density 
perturbations and dissipative forces for distant 
inclusions 

In what follows we consider a non-equilibrium steady- 
state problem by setting d T n = 0 in Eq. 

V 2 n - UV 2 n + nX7 2 U - (g, V)n(l - n - U) = 0, (11) 

where obstacles are given by a distribution U of the heavy 
gas component. Far from small obstacle (whose size is 
comparable with lattice constant) the density distribu¬ 
tion n = no + Sn weakly deviates from the equilibrium 
one no [45]. In this case, interaction between gas parti¬ 
cles is less significant and the drift term in Eq. © can 
be written in the linear approximation n 2 ~ rip + 2noSn. 

Simple analytical expressions for density perturbation 
and dissipative forces for the ensemble of widely sepa¬ 
rated small obstacles can be obtained using the qualita¬ 
tive approach described in Appendix [A] This approach 
is similar to that based on the method of molecular field 
that was used to describe the elastic interaction of col¬ 
loidal particles in a liquid crystal, see jHg]. In particular, 
the gas density perturbation far from an isolated obstacle 
can be written as 

«5n(r)~(n,V P )G(r), ( 12 ) 

where G(r) is anisotropic screened Coulomb-likc poten¬ 
tial that takes the form 

1 p - 9 |r|+qr 

CM = (13) 

in 3D case, and 

G(r) = ( 14 ) 

in 2D case. Here, /\ 0 is the modified Bessel function, 

vector q = (1/2 — n o)g determines the preferable direc¬ 
tion of screening and depends on external sweeping field 
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g (or gas flow) and on equilibrium gas concentration no 
(bath fraction), ft plays a role of the molecular field or 
an average flux near obstacle, see Appendix lAl 

At low concentrations of gas (no < 1/2), the dense re¬ 
gion of the gas ahead of the inclusion is described by an 
exponential asymptotics, while the asymptotics of the de¬ 
pletion region behind the inclusion is power-law. When 
gas concentration increases and no becomes greater than 
1/2, the anisotropy vector q = (1/2 — no)g changes its 
direction to the opposite. It means that switching of 
the wake direction occurs together with corresponding 
switching between the exponential and power-law asymp¬ 
totics. The distribution Sn( r), related to the anisotropic 
screened Coulomb-like potential, formally describes a 
“medium polarization” around the inclusion induced by 
an asymmetrical “dipole” (see, e.g., Fig. [TJ. 

In the case of ensemble of distant inclusions, the force 
exerted by the ith inclusion on the fctlr one can be roughly 
estimated as (see Appendix |A} 

fki ~ -V Rfc (Hi, V R J G(R fc - Ri), (15) 

where R*, is the center of the A;th inclusion. The 
anisotropic screened Coulomb potential G, giving the 
asymmetrical form of the obstacle wake m, naturally 
leads to the non-Newtonian character of dissipative forces 
acting between obstacles, fki —fik- As seen from 
Eqs. m and csd, the asymptotic behaviors of the den¬ 
sity perturbation and the dissipative forces acting be¬ 
tween widely spaced small inclusions are defined by the 
moments of a screened anisotropic Coulomb potential. 
The local density perturbation around an obstacle is 
formed by an effective flow 17 (molecular field) that is 
determined by the external flow and flows induced by 
gas density perturbations of all the inclusions. The latter 
means that the interaction between two inclusions cannot 
be separated out of the influence of all other inclusions. 
This is a general property of a nonlinear response or non¬ 
linear systems, see [581] . The employed approach allows 
one to consider the non-linear response and to represent 
expressions for Sn and fkj in the form similar to that 
given by the linear response for moving probe-particles 
(cf. note [HI), the only difference is that asymptotic be¬ 
haviors are associated with anisotropic screened Coulomb 
potential instead of Coulomb one, |r| -1 , and with mean 
gas flow (mean field) near inclusion instead of velocity 
of a probe-particle. However, Eqs. m and m are ob¬ 
tained within somewhat naive approach and give only a 
qualitative picture of dissipative interaction. 

More rigorous results for the asymptotics behavior can 
be obtained within the single-layer potential approach 
for inclusions with sharp boundaries. Representation of 
solution for Sn in the form of single-layer potential was 
proposed in [45( to describe the gas density perturba¬ 
tion around single obstacle in 2D case. In this paper, 
we use this representation and its multipole expansion 
(Appendix [B]) to find a general form of asymptotic be¬ 
havior of dissipative forces for widely separated obstacles. 
Particularly, in 3D case this method gives the following 


asymptotic behaviors: 

g-fld+qr _ 

5n(r) -—-J(r,q) (16) 

for density perturbation caused by a small isolated ob¬ 
stacle and 

e -g|rfci|+qr fci r 

fki ~~~—rn-i— / n (xfc)/(r fe j,q,x fe ) dx fe (17) 

47T | r fci | J Sk 

for the dissipative force exerted by ith inclusion on 
kth one in the dipole approximation, when the distance 
|rj/-| = |Rj — Rfc| between inclusions is much larger than 
their radii a^a/j) ~ £, n(x^) is the exterior normal at 
the point x& on the surface of the fcth inclusion. For 
simplicity, we have considered spherical obstacles, S/ is 
the surface of the £:th inclusion (|xfc| = etfc). Functions 
/(r, q) and I (r ki, q, x^) have a power-law dependence on 
l/r and 1/ru, respectively [see Eqs. (IB26I) and (IB21I) ]. 
In particular, function / (r ki, q, x*,) can be represented in 
general form as 


I = A(q,a l ,x fc ,) + B(q,o i ,x fc ,) (q -q 


r ki 

I Dei I 


rfci 

|rfeip 


(18) 

where A and B are determined only by the obstacle sur¬ 
face and external field g. 

In 2D case, Sn and fki are determined by the potential 
exp(qr)A' 0 (gr), Eq. (THl) . having the asymptotic behavior 
~ r -1 / 2 exp(qr — qr ) at large r. Detailed form of the 
density distribution Sn around a single circular obstacle 
have been considered in [45J. The leading asymptotics of 
the dissipative force and its comparison with numerical 
results for Eq. ([!]) in 2D are given in Appendix iBl 

In the particular case of half filling (no = 1/2), q = 0 
and the potential G degenerates into usual Coulomb one, 
see [45[ and Appendix [B] The form of the interaction be¬ 
tween obstacles corresponds to anti-Newtonian dipole- 
dipole one as it is in the case of the linear response 
j28| . Note that single-layer potential approach enables 
not only correct description of the dipole-dipole interac¬ 
tion but also accounting for the higher order multipole 
moments. 

The linear flow approximation allows us to describe 
the concentration-dependent switching of dissipative in¬ 
teraction between obstacles and to determine the type 
of this interaction. The latter, being formally expressed 
by Eqs. m or 03, belongs to the induced dipole- 
dipole (generally multipole) type of interaction in the 
non-equilibrium steady-state case. In contrast to usual 
electrostatic interaction between polarizable particles in 
electric field, Eqs. m and G3 describe the interaction 
between induced “asymmetric” dipoles (with nonzero to¬ 
tal induced “charge”, see Appendix [Bit, that is associated 
with anisotropic screened Coulomb-like potential with 
preferential direction of anisotropy q. This approxima¬ 
tion is valid for small and widely separated obstacles and 
does not describe the non-linear effects that are essential 
for closely located obstacles as well as in the vicinity of 
a large-sized obstacle surfaces. 
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Figure 2. Stationary wakes (numerical results) with kink¬ 
like profiles Sn(x, y = 0) that describe cavities (a) ahead of 
and (b) behind the obstacle; no = 0.3 for (a) and no = 0.7 
for (b). (c) Concentration profiles n(x,y = 0;no) at sev¬ 

eral values of the equilibrium concentration no- Note that 
a compact jammed region grows ahead of the obstacle as no 
tends to 1/2, for no exceeding the half filling the wake profile 
becomes inverted, (d) Contour plot of concentration distri¬ 
bution 5n(x, y) = n(x, y) — no for no = 0.38. Here, |g| = 0.5, 
u = 1, a = 7 (in units of f), and spatial coordinates are in 
units of i. 


B. Nonlinear blockade effect near surface of big 
obstacle 

Here, we briefly discuss the non-linear effects caused 
by the gas particles blockade resorting to the numeri¬ 
cal stationary solution of two-dimensional equation © 
in the mean-field approximation. For a relatively large 
obstacle, whose size is much larger than the lattice con¬ 
stant, the screening of the gas flow near the obstacle sur¬ 
face leads to a growth of the obstacle’s effective size. As 
a result, a compact high-density (jellium-like) region is 
formed ahead of the obstacle, Figs. ©a), (c), and (d). 
Figure ©a) shows that behavior of 5n near the obstacle 
surface has a pronounced step-like character at no < 1/2. 
Note that such a step-like behavior of the density per¬ 
turbation Sn is quite expected since the general type of 
equation © admits kink-like solutions, e.g., for a one- 
component lattice gas, u = 0. As gas concentration no 



Figure 3. Kink-like concentration profile n(x, y) formed ahead 
of two closely located obstacles at no = 0.3. The distance be¬ 
tween the obstacle centers equals 4a, other system parameters 
are the same as in Fig. [5] 


approaches 1/2, the compact dense region grows (while 
its boundary becomes diffused), see Fig. ©c), until the 
uniformly decreasing distribution is formed. Further, at 
no > 1/2, the upstream part of the profile transforms 
into an inverted diffusive wake with an extended dense 
region ahead of the obstacle, Fig. ©c), while a local¬ 
ized low-density region (resembling the form of a cavity) 
with an inverted step-like profile is formed downstream, 
Fig. ©b). Note that a similar compact structure occurs 
in a dusty plasma f© . Il3l ] . That structure is formed by 
a flow of smaller dust grains ahead of a void formed by 
larger grains. 

A similar nonlinear effect occurs for closely located in¬ 
clusions when their individual density perturbation coats 
are overlapped. The overlapping leads to an additional 
screening of the gas flow and to the formation of a com¬ 
mon non-linear coat around them with a step-like behav¬ 
ior of the density perturbation profile, Fig. © at least 
in 2D case. Note that formation of a common coat can 
signify the effective pairing between the inclusions (at 
no > 1/2), i.e., formation of a stable coupled doublet 
[4J] (see also subsec. IIII Cl) . 


C. Dissipative force switching for two moderately 
separated big obstacles 


We next consider numerically the wake-mediated force 
between two obstacles for two orientations of the line 
of their centers — parallel and perpendicular to the gas 
flow. We use Eq. © in the mean-field approximation 
that takes into account the non-linear blockade effect for 
gas particles. The total force exerted on a given obstacle 
includes the part associated with the individual friction 
force and the one associated with the influence of another 
obstacle. To separate out the inter-obstacle contribution 
from the total dissipative force we consider the quantity 
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Equilibrium concentration n 0 (fractions of 1) 


Figure 4. Transverse alignment. Dependencies of y- 
components of dissipative forces /12 and /21 against the equi¬ 
librium gas concentration no- Several regimes of the drive 
|g| = 0.3; 0.4; 0.5 are plotted for comparison. Other system 
parameters are the same as in Fig. [l] the forces are in units 
of kT/£ (£ is the lattice constant). 



fij = fi - ft = f [Mr, R-i, Rj) - Sn( r, R*)] Vuj(r) dr, 

(19) 

where fi is the total force acting on the ith obstacle in the 
presence of the jth one and /f is its individual friction 
force. 

Transverse alignment (Fig. [I]). From the symmetry 
of this configuration it follows that the y-components of 
the forces two obstacles exert on each other are equal 
and opposite, f\ 2 = — f'f ■ At low equilibrium concen¬ 
trations (no < 1/2), Fig. Ufa), the dissipative interaction 
manifests itself as an effective repulsion between the ob¬ 
stacles, since f 2 1 < 0 and ff 2 > 0, see Fig. [U Qualita¬ 
tively, this effective repulsion is simply explained by the 
overlap of the density coats around the obstacles that 
leads to formation of a dense region between them act¬ 
ing like a repulsive barrier, see Fig. [Ha). In contrast, at 
no > 1 / 2 , the overlap of the individual density pertur¬ 
bation coats of the obstacles results in formation of an 
extended dense zone ahead of them that blocks the gas 
flow, so that the region between the obstacles becomes 
depleted. As Fig. [4] suggests, this collective blockade ef¬ 
fect of gas particles leads to effective attraction between 
obstacles in a dense medium, f 21 > 0 and ff 2 < 0 . 

Thus, when gas concentration n 0 increases, the dissi¬ 
pative interaction between the obstacles switches from 
effective repulsion to attraction. In the no = 1/2 case, 
the effective interaction between the inclusions vanishes, 
fi 2 = f21 = 0 , regardless of the distance between them. 
The dissipative interaction between the inclusions natu¬ 
rally vanishes in the limit of empty medium —> 0, due 
to wake depletion. The same is true in the total jamming 
limit ?iq —> 1 . 


0 n 0 1 



-50 0 50 100 150 



-100 -50 0 50 100 



-150 -100 -50 0 50 


Figure 5. Longitudinal alignment. Steady-state concentra¬ 
tion distributions (average occupation numbers) n(x, y) of the 
gas particles near the obstacles, evaluated numerically within 
the mean-field approximation of Eq. Q. corresponding to 
the equilibrium concentrations no = 0.2 (a), 0.5 (b), and 
0.8 (c). Panels (a) to (c) illustrate three different regimes of 
the dissipative interaction, see Fig. [6] (a) effective attraction 
(I/21I > I/ 12 I), (b) anti-Newtonian interaction (f 21 = /f 2 ), 
(c) effective repulsion (|/fj| < |/i 2 |). The external field g 
(|g| = 0.5) is directed along the z-axis; the impermeable 
(u = 1) circular obstacles are of radius a = 7 (in units of 
£), their positions are marked with the black circles; the dis¬ 
tance between the obstacles’ centers equals 10a. The gray 
background corresponds to the equilibrium gas concentration 
no for every contour plot, in consistence with the colorbar; 
spatial coordinates are in units of £. 


Longitudinal alignment (Fig. [5]). At low concentrations 
(no < 1/2), a typical situation for Brownian systems 
takes place: An inclusion falling on the depleted wake 
induced by another inclusion is effectively attracted to 
it since the friction force in depleted regions is weaker 
0 H3- This type of effective interaction is often re¬ 
ferred to as the wake-mediated [ 59 J, As Fig. [5](b) 
suggests, the second obstacle does not practically affect 
the first one, /f 2 ~ 0. In contrast, at high concentrations 
(no > 1 / 2 ), the second obstacle does not feel the influence 
of the first one, f 21 ~ 0 , whereas the first obstacle comes 
under the excess pressure of the dense gas region created 
ahead of the second one due to the blockade effect. As a 
result, the effective interaction changes its sign, switching 
from effective attraction to repulsion, Fig. (Hl[b). In the 
case of no = 1 / 2 , the effective interaction between the in¬ 
clusions becomes strictly anti-Newtonian, /f 2 = /fj 7 ^ 0, 
see Fig. [U[b). Note that for a dense gas in the blockade 
regime, the second obstacle “pushes” the first one up¬ 
stream, thus reducing the total friction force /f exerted 
on the first obstacle, Fig. [HI a). 
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0.4 



Equilibrium concentration n 0 (fractions of 1) 


Figure 6. Longitudinal alignment. Concentration dependence 
of the total forces /f (no) and /if (no) acting on each obstacle 
(a), and the forces /i 2 (no) and /li(no), acting between the 
obstacles (b), at three magnitudes |g| = 0.3, 0.4, 0.5 of the 
drive field. Other system parameters are the same as in Fig. [5] 
forces are in units of kT/l (£ is the lattice constant). 


The above described behavior of forces, Fig. [Gl can be 
qualitatively explained by using the results of the lin¬ 
ear flow approximation. For example, for point-like ob¬ 
stacles at no < 1/2, forces /f 2 and /f 1; see Eq. (fl5l) . 
are associated with potentials oc exp(— 2 qr\ 2 )/^/ri 2 and 
oc 1/y / ? 2 i, respectively [see asymptotic expression (IA15I) 
in Appendix IX], so that |/f 2 | <C |/fil- Besides, the single¬ 
layer potential method gives correct leading asymptotics 
f\i ~ |ri 2 1 -3 / 2 at large |ri 2 |, that is in satisfactory agree¬ 
ment with numerical result for the general non-linear 
problem, Eq. m , see Appendix iBl 

Note that for closely located obstacles the non-linear 
inter-obstacle attraction can determine the dissipative 
pairing by the creation of common perturbation coat 
around them. The effect of a similar nature was obtained 
earlier in [24| for two driven tracers. Indeed, at high gas 
concentration (no > 1/2) the depleted cavities formed 
around each obstacle, see Figs. 0c) or 0c), can entail 
specific behavior of dissipative forces depending on the 
distance between obstacles. In particular, the effective 
interaction between two obstacles in close proximity un¬ 
dergoes an abrupt change in the asymptotic behavior, see 
Appendix [B] and figures therein, that can be indicative 
of the dissipative pairing effect. 

Dependence of the strength of dissipative interaction 
between obstacles on the external driving field g appears 



Figure 7. Dissipative forces (a) /fi (for longitudinal align¬ 
ment) and (b) / 21 (for transverse one) versus the external 
drive g for the gas concentrations no = 0.1; 0.2; 0.3. Obsta¬ 
cle sizes and separation correspond to that on Figs. [T] and [5] 
Schematic illustrations of the shape transformations of den¬ 
sity perturbation coats are shown for transverse alignment. 

to be non-monotonic, see Fig. [7] The characteristic peak 
of interaction corresponds to the drive magnitude when 
the most efficient common density coat is formed around 
the obstacle pair. This behavior can be explained by the 
changes in the shape of density perturbations, e.g., for 
the case 0 < no < 1/2, see Fig. [3(a). At low gas concen¬ 
tration the effective repulsion between obstacles vanishes 
in the limit of weak driving, since slow flow of a sparse gas 
does not induce significant gas perturbations and, thus, 
wake-mediated interactions. The characteristic peak of 
interaction corresponds to the driving magnitude when 
the common density coat is formed (see schematic illus¬ 
trations on Fig. 0: in this regime, profile of the den¬ 
sity perturbation provides the most efficient dissipative 
wake-mediated influence between obstacles. Strong driv¬ 
ing field causes the perturbation coat around each ob¬ 
stacle to decrease in lateral dimension and increase in 
longitudinal, so that overall density coat stretches along 
the flow direction. As a result, overlap of the individual 
obstacles’ coats reduces, and their mutual influence de¬ 
cays. In other words, strong enough drift flow reduces 
the common density coat. This qualitative reasoning is 
also true in the case of effective attraction under longitu¬ 
dinal alignment. Note that the peak position shifts and 
increases towards the region of strong driving as gas con- 
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centration decreases in case of longitudinal alignment of 
obstacles, Fig. [7Jb), while in case of transverse one the 
situation is just the opposite, Fig. [7|a). Hence, the most 
favorable condition for the pronounced common coat or¬ 
ganization is determined by both the equilibrium gas con¬ 
centration no and the strength of external driving field 
g- 

The magnitude of the evaluated forces can be easily 
estimated, e.g., for the case of atoms adsorbed on solid 
surface. Choosing the lattice spacing parameter to be 
£ = 3 A, at room temperature one obtains the range 
of dissipative forces to be 5-10 pN [see Fig. [Htb)] , while 
the friction force is approximately one order of magni¬ 
tude stronger [see Fig. [Oja)]. Notice that the same ratio 
between magnitudes of friction and dissipative forces is 
observed for probe-colloids moving through a colloidal 
suspension in 3D case [1, |b|. In addition, as is seen from 
Fig. [3a), at concentrations close to no = 1/2 the forces 
exerted on each obstacle by the gas are almost equal, i.e., 
the dissipative interaction between obstacles takes anti- 
Newtonian character, /f 2 = /|j, see Fig. [3b). It should 
be mentioned that an analogous behavior occurs for two 
probes moving along their line of centers through a col¬ 
loidal suspension: both inclusions may experience the 
same drag force, as was observed in a recent experiment 
Q, at the effective volume fraction of 0.41. However, in 
this case, the effect is due to the hydrodynamic interac¬ 
tions between bath particles. 


IV. CONCLUSION 

Let us briefly summarize the main obtained results for 
dissipative (wake-mediated) interaction between obsta¬ 
cles embedded into gas flow with taking into account 
short-range repulsive interaction between gas particles. 

- We have shown that increasing of the gas concen¬ 
tration enhances the role of inter-particle interaction and 
can lead to the sign change of the effective interaction be¬ 
tween obstacles, i.e., switching from effective attraction 
to repulsion or vice-versa. The effect of concentration- 
dependent force switching is associated with obstacle 
wake transformation — its inversion. 

— In the case of small and widely separated obstacles, 
the wake-mediated dissipative interaction between them 
has been shown to belong to the type of induced dipole- 
dipole (generally, multipole) interaction associated with 
anisotropic screened Coulomb potential. To this end, we 
have developed the representation for the gas density per¬ 
turbation in the form of single-layer potential. Formally, 
this is a generalization of the single-layer potential ap¬ 
proach for the electrostatic interaction between polariz¬ 
able particles induced by stationary external field. Our 
approach is applicable to non-equilibrium steady-state 
case where interaction between the obstacles is induced 
by gas flow. Obtained analytical expressions qualita¬ 
tively explain the asymmetry of the obstacle’s wake, the 
long-range behavior of dissipative interaction, its non- 


Newtonian character, and switching of both the wake 
direction and the dissipative forces. 

- Dissipative interaction between obstacles is most 
pronounced when a common perturbation coat around 
them (collective wake) is formed. The force depends non- 
monotonically on equilibrium gas concentration, magni¬ 
tude of external sweeping field (gas flow), and alignment 
of the obstacles. In particular, at low gas concentrations 
two obstacles are effectively attracted in the case of lon¬ 
gitudinal alignment and repel each other in the case of 
transverse one. At high gas concentrations the situation 
is just the opposite. 

— The non-linear blockade effect of gas particles is sig¬ 
nificant near the surface of relatively big obstacles and/or 
for closely located ones. In this case, repulsive interaction 
between gas particles has been shown to lead to screening 
of the gas flow near the obstacles and to formation of a 
common coat of gas density perturbation around them, 
with pronounced step-like behavior of the density profile. 
Formation of common coat can determine the non-linear 
mechanism of dissipative pairing between the obstacles 
(see, e.g., @13). 

It should be noted that we initially used rough ap¬ 
proximations, so that a number of important questions 
were left behind the scope of our paper. In particular, 
using the mean-field approximation, we lose information 
on the short-range correlations in a gas, such as “back 
correlations”, see, e.g., [36|,[3; which have to occur near 
the obstacle surfaces (another gas component). In addi¬ 
tion, neglecting fluctuations in a gas, i.e., the term (5J“ 
in Eq. ©, we do not take into account the fluctuation- 
induced (Casimir-like) forces, see, e.g., [H, l6ll - t66j which 
can be significant for pairing effect at small inter-obstacle 
distance. 

Obtained results may be of interest when considering 
the dissipative structure formation (see, e.g., ®) , collec¬ 
tive friction force or collective energy losses in an ensem¬ 
ble of inclusions, and can find applications in systems 
with driven hop ping transport (e.g., surface kinetics of 
adsorbed atoms [U [23|, [3 jbj], fast ionic conductors, 
etc.) or serve as a rough model for colloidal suspensions 
or dusty/complex plasma EMI- 
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Appendix A: Qualitative picture of dissipative 
interaction: a rough analytical approach 

In this appendix, we roughly estimate the wake- 
mediated interaction between widely separated small ob¬ 
stacles imbedded into gas flow. Let us consider a non¬ 
equilibrium steady-state problem in the long-wavelength 
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approximation, Eq. m : 

V 2 n - t/V 2 n + nS7 2 U - (g, V)n(l - n - C7) = 0, (Al) 

where inclusions are given by a distribution U of the 
heavy gas component. For simplicity, we consider a 
smooth distribution U(r) = u (r — R fc ), where distri¬ 
bution u(r — Rfc) describes fcth inclusion and has a com¬ 
pact carrier located near the inclusion center Rfc. For dis¬ 
tant inclusions, we assume that f u( r—Rfc)tt(r -Rj)dr w 
0 . 

For widely separated small obstacles (whose sizes are 
comparable with the lattice constant), interaction be¬ 
tween the particles is less significant and the drift term 
in Eq. ED can be written in the linear approximation, 
see [45[ . Assuming that distribution n = no + Sn weakly 
deviates from the equilibrium one no, we linearize the 
drift flow term in Eq. ED, taking n 2 « n§ + 2noSn,, and 
rewrite the equation in the following form: 

V 2 Sn — 2(q, V)Sn 

= UV 2 6n - (n 0 + Sn)V 2 U - (g, V)(n 0 + Sn)U, (A2) 

where q = (1/2 — no)g. Based on Eq. (IA2I) . we estimate 
the asymptotic behavior of the dissipative interaction be¬ 
tween the obstacles depending on the distance between 
them, their mutual alignment, and equilibrium gas con¬ 
centration no- We shall use a qualitative approach which 
allows us to obtain simple analytical expressions for den¬ 
sity perturbation and dissipative forces. 

It is convenient to consider an integral representation 
of Eq. (IA1I) using the Green’s function G(r — r') of the 
equation 

V 2 G(r - r') - 2(q, V r )G(r - r') = -6(r - r'). (A3) 

The form of this Green’s function is similar to the 
anisotropic screened Coulomb potential 

p —g|r—r'|+q(r—r') 

G(r - r') =-—--- (A4) 

47r|r — r'| 

in 3D case and 

G(r - r') = ^-e q(r ~ r ) AT 0 (g|r - r'|) (A5) 

Z7T 

in 2D case. By using (1A3ll . we can rewrite the equation 
for the gas density perturbation Sn in the form 

Sn( r) = [no + 6n(r)]?7(r) 

+ J U(r') (fi(r'), V r ) G(r — r') dr', (A6) 

where 

fl(r) = 2g(l - no)[no + Sn( r)] - 2V r <5n(r). (A7) 

Equation (IA6I) can be simplified by applying an approach 
similar to the self-consistent molecular field approach 
[58^ . Since distribution U of the heavy component is 


localized near the inclusion centers Rj and has compact 
carriers Uj( r) = u( r R j ) < 1, we may consider Uj( r) 
as a probability de nsity dist ribution and the integral in 
ESI) as an average (H, VG) ■ associated with this distri¬ 
bution. Here, (...) ■ = /(.. .)u(r — Rj) dr. Then, using 
the mean-field approximation, (fl, VG) • ~ (flj, VGj) « 
(Dp VG(r - R ; )), Eq. EH) can be rewritten as 

Sn( r) r: [no + Sn(r)]U(r) + ^ ^fl(Rj), V r ^ G(r - Rj), 

(AS) 

where 

fi(Rj) = 2g(l — no) [no + Snj\ — 2 X78nj (A9) 

plays the role of a molecular field or an average flux in the 
system, these quantities being defined by external field 
g and the density perturbation field due to other inclu¬ 
sions. Equations for the constants Snj and S75nj can be 
obtained in a self-consistent manner by using Eq. (1A8[) , 
see Ref. [68j . 

Representation ESP enables us to estimate qualita¬ 
tively the asymptotic behavior of gas density perturba¬ 
tion far from an isolated inclusion and the asymptotic 
behavior of the dissipative force between widely sepa¬ 
rated inclusions. Using (IA4I) . gas density perturbation 
(IA8I) far from an isolated inclusion can be written as 

Sn(r) ~ (n,V P )G(r), (A10) 

that is 

_ p -«M+qr 

fc(r)~ (ft,V r )-—- (All) 

in 3D case and 

_ p -«M+qr 

Mr)~(^V r )^i^- (A12) 

in 2D case. In the last expression, the asymptotic be¬ 
havior of the Bessel function /<o(< 7 |r|) ~ |r|~ 1 / 2 e -9 l r l 
for large r was used. At low concentrations of gas 
(no < 1/2), the dense region ahead of the inclusion 
is described by an exponential asymptotics, while the 
asymptotics of the depleted region behind the inclusion is 
power-law. When gas concentration increases and no be¬ 
comes greater than 1/2, vector q = (1/2 — no)g changes 
its direction. It means that switching of the wake direc¬ 
tion occurs, together with corresponding switching be¬ 
tween the exponential and power-law asymptotics. At 
no = 1/2 we have q = 0 and the asymptotics of pertur¬ 
bation Sn corresponds to a dipole-like polarization of gas 
density perturbation around the inclusion. 

The force acting on the fcth inclusion is fk = 
J <5n(r)V r u(r — Rfc) dr ~ —\75nk- For small (point-like) 
inclusions, the force exerted by the jth inclusion on the 
kth one can be roughly estimated as 

h, - -V R , (n,-,V Rfc )G(r), 


(A13) 
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that is 

_ e -g|Rfe-Rj|+q(R-fe-Rj) 

fkj -V Rfc (Slj, V Rfc )- | Rfc _ R | - (A14) 

in 3D case and 

_ e -9|R-fe-R-3l+q(R-fe-R-i) 

fkj ~ -V Rfc (n jt V Rfc ) | Rfc _ R .j_ 1/2 . (A15) 

in 2D. As is seen from Eqs. (IA10I) and (IA13I) . the local 
density perturbation around an obstacle is formed by an 
effective flow ft (molecular field) that is determined by 
the external flow and the flows induced by gas density 
perturbations of all the inclusions. 

Appendix B: Single-layer potential approach 


where Q n = (Q,n r ), q n = (q, n r ), n r is the outward 
normal at the point r € 5), e + = 1 outside the inclusions 
and £~ = (1 — Ui) 2 inside the ith inclusion, and notation 

Vn (• • ■) = li m |F— Ri |—Hr—Ri|±o (%r) (?) is used. 

The solution of (IB3I) (IB4I) can be represented in the 
form of a single-layer potential, similarly to that used in 
[4f| for a single obstacle, 

Sip(r) = '^2f G( r-r , )At»(r , ) dr ') ( B5 ) 

i Js * 

where G(r — r') is the Green’s function, (IA4I) in 3D, 
or (IA5I) in 2D. The quantity m{ r') plays the role of a 
“charge” density induced by the external field g on the 
obstacle surface Si (6{|. It satisfies the following integral 
equation determined by the matching conditions (IB4I) : 


A more rigorous result for wakes and for the dissipative 
force can be obtained in the framework of the single-layer 
potential method for inclusions with sharp boundaries. 
Again, we start from Eq. (fill) . It is convenient to use 
a new function D(r) = n(r)/[l — -u(r)] governed by the 
equation (see [45f l 

V {e [V^ — V 7 (1 — ip) g]} = 0, (Bl) 

where e = e(r) = [1 — it(r)] 2 , and it is assumed that 
u( r) / 1. Let us represent the solution tp(r) m i/j 0 + 
5ip( r) as a small deviation 5ip(r) from the equilibrium 
distribution ipo = ng, and linearize Eq. EH): 

V [e (Vch/> - 2q Sip - Q)] = 0, (B2) 

where Q = no(l — no)g and q = (1/2 —n 0 )g. This linear 
equation takes into account the interaction between gas 
particles in the first order of the perturbation theory. In 
this sense, (IB2I) is the simplest possible generalization of 
the drift-diffusion equation that was exploited in [27t for 
a gas of non-interacting particles at low concentrations. 

The inclusions are represented by the distributions of 
heavy gas-component Uj( r) = u( r — R,) centered at 
points Rj with homogeneous concentration Uj( r) = Uj = 
const inside inclusions and Uj (r) = 0 outside them. Note 
that in the case of inclusions with sharp boundaries, 
Eq. m allows for a solution in the class of continu¬ 
ous functions, whereas function n(r), obeying Eq. ED, 
as well as its normal derivative have a jump at the inclu¬ 
sion’s boundary. The density perturbation inside (Sip-) 
and outside (5ip + = Sn) the inclusions obey the equation 

V (VdV’ ± - 2q5ip ± - Q) = 0. (B3) 

Equation (IB3I) is supplemented by the matching condi¬ 
tions for 8ip± on the surface S t of ith inclusion: 

5ip + ( r) = 5ip~( r), 

e + [V+5V + ( r ) - 2q n 6i/j + (r) - Q n ] 

= e T [V'<5^ _ (r) - 2q n 5il>- (r) - Q n \ , (B4) 


2A, [V+ - 2q n (r t )] ^ 

3 


Clr, i - , 1 dr, 


+(A* — l)ni(vi) = 2XiQ n (ri), (B6) 


where r* e Si and A, = A (Hi) = (e + - e~)/(e + + e“). 
Equation (1BGI) was derived with the use of the jump 
theorem for the normal derivative of the potential of a 
single layer on an obstacle surface [7Cj, V„(5'i/> ± (r) = 
=F/i(r)/2 + V n (5^>(r). Representation (IB5I) and Eq. (IBGI) 
describe the general solution for obstacles with arbitrary 
geometry of their surfaces (Lyapunov surface, see 0 )- 
Considering that r) = Sn( r) and using Eq. (flOl) . 
we can write final expression for the density perturbation 
around obstacles and the force acting on an obstacle, 
both being induced by sweeping field g (or by the gas 
flow): 


Sn{r) = 

3 


G(r — rj )fj,j (rj ) div,-, 


(B7) 


fk = ~ Y / n ( r fc)G(rfe — drfcdr^. (B8) 

j Js k J Si 

This representation of the solution has direct analogy 
with induced interaction between dielectric particles in a 
stationary electric field Q: External electric field induces 
charge n on the particle surface, leading to its polariza¬ 
tion, e.g., inducing the dipole moment for a spherical 
particle, see Q. This, in turn, leads to multipole (e.g., 
dipole-dipole) interaction between the particles. How¬ 
ever, in our case, contrary to the electrostatic problem 
the density /rfc is induced by an external field (flow) 
on the inclusion surfaces, and multipole interaction be¬ 
tween them is determined not by the Coulomb poten¬ 
tial |r| _1 but anisotropic screened Coulomb-like poten¬ 
tial |r| _1 exp(qr — g|r|) (in 3D case), see (IA4D and (IA5I) . 
Such form of the potential leads, in particular, to noncon¬ 
servation of induced surface density, J /idS 0, and to 
asymmetric distribution of “induced potential” Sn near 
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the inclusion. The latter describes inclusion wake, e.g., 
wake with a localized region of dense gas ahead of the in¬ 
clusion and an extended depleted tail [see, e.g., Fig.QJa)]. 

In the particular case of the half filling (no = 1/2), 
the second term in equation (IB2I) vanishes (q = 0) and 
the problem is reduced to an electrostatic-like problem 
V [e(V(5 ip — Q)] = 0 for dielectric particles in a uniform 
electric field Q = g/4. In this case, density distribution 
fo/>(r) is similar to the distribution of the electrostatic 
potential characterizing the scattered field. It means 
that the induced interaction between obstacles via their 
common environment (density perturbation) behaves like 
electrostatic dipole-dipole (generally, multipole) interac¬ 
tion. For a single obstacle with radius a, density pertur¬ 
bation Sn = 5 'i/j + around the obstacle at no = 1/2 can be 
obtained in an explicit form: Sn = Aa 2 (Q, V r ) In a|r|“ 1 
for 2D case and Sn = Aa 2 (Q, V r )|r| _1 for 3D. These 
results explain both the power-law asymptotic behav¬ 
ior of gas perturbation and the symmetry of the “up¬ 
stream/downstream” tail, see Fig. EKb) and [45|. This 


case (no = 1/2) corresponds_to the linear response of Sn 
to the external field g, cf. (28|. Note that symmetry of 
wake (or profile of perturbation) generated in a medium 
by a moving probe particle is a common result for sys¬ 
tems described in the linear response approximation (see, 

e.g., (3). 

For widely separated inclusions, when the distance be¬ 
tween their centers | r^^ | = |Rfc — Rj| is much larger than 
their characteristic sizes dj, \r k j\ = |Rfc — Rj| a k (dj), 
the multipole expansion of the potential G can be used: 

G(r k - rj) « G( r kj ) + (x kj ■ V rfcj .) G(r kj ) + ■ ■ ■ , (B9) 

where x k j = x k — x, and x k = r k — R&. Next we 
consider the particular 3D case for spherical obstacles 
with radii a k . For obstacles located far from each other, 
|rfcj| |xfcj|, one can use the multipole expansion (IB9I) 
for the kernel of integral operator in (IB6I) . In the dipole 
approximation, the integral equation for the induced sur¬ 
face density fj, k (x k ) on the surface of the kth inclusion 
takes the form 


E e -q\r k i\+Vki r 

——-:- / (1 + XfcjUfcj) ptj(Xj) dx,; = x fc Q, 

i^ k ^ n \ r ki\ J Si 


where 


Ufci(r fc j) = q -q 


^ki I *ki 


rfci r ki 


and A Xfc denotes the integral operator for a single obstacle 

A Xk Mfc(xfc) = a k '— — t- tfc(xfc) + x fc (Vj fc - 2q) / 

J s 


3 -g|xfcfc|+qx fcfc 


Sk 4w\x kk \ 


■ Mfe(xfc) dx' fc . 


(BIO) 


(BH) 


(B12) 


Equation (IB 101) for jx k has small parameter exp (— q\r ki \ + qrfcj)/47r|rfcj| -C 1, that allows us to consider the influence 
of other obstacles on a given one as a small perturbation fi l k of the solution fi° k = pi° for a single obstacle, /j, k ~ m/ + a 4- 
In this approximation, equations for = pt° and (i\ take the form 


A Xfc Mfe(xfc) — XfcQ, 


(B13) 


AMfe(xfc) = -x fe (V Xk - 2q) ■ 


i^k 


-Q’kjMl+qrM 

47r|r fci | 


I Si 


(1 + x fci u fei )/i°(x.j)dx i . 


The formal solution of the last equation can be written as 


Mfc( x fc) = J2 


o~<lWki\+Vki 


i^k 


47r|r ki 


-Ax, 1 / W (xfc, Xi, r k i) pt°(xj) dxi, 


where 


W'(?C-k X-i i Yki) — 2(c[X/ c )(l Xfc^ U./^ ) 


(B14) 


(B15) 


(B16) 


The dissipative force acting on an obstacle is deter- mined by the gas density perturbation Sn( r) on its sur- 
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face, Eq. m- In this case we can set Sn ss Sip(r ), where 
r » R fc + Xfc and |xfc| = a k is the radius of the fcth ob¬ 
stacle. In the dipole approximation, Eqs. (1B9I) . (IB 131) . 
(IB14I) , density perturbation near the fcth obstacle can be 
written in the form 

Sn k « Sn° k + [ G(x fcfe )/4( x fc) dx' fe 

Js k 

+ E f [G'( r fci) + (xfci-V rfci )G(r fci )]M°(x' fc ) dx' fe . 

ijik JSi 

(B17) 

The right-hand side of expression (IB17I) containing the 
sum over all the obstacles i ^ k describes their direct 
influence on the given fcth obstacle. The first term in 
Eq. (1B17I) 


Sn°k = [ G(x fcfc )/ig(xfc) dx' fc (B18) 

Js k 

gives the contribution to the gas perturbation around the 
fcth obstacle caused by the fcth obstacle itself. 

Using Eqs. (IB 151) and (IB17D . contribution to the den¬ 
sity perturbation Sn k near the fcth inclusion caused by 
other inclusions in 3D case can be written as 


Sn k 


Sn °k ~ E^’ 

i^k 


where 


5n k i 


e ~q\*ki\+<Vki 

47r|r fei | 


I (rfcj, q, Xfc) 


(B19) 


(B20) 


is contribution of the ith obstacle to the density pertur¬ 
bation near the fcth obstacle surface, 


Expressions (IB20I) (IB22I) are obtained in the dipole ap¬ 
proximation and give a rough asymptotic behavior of 
the induced non-equilibrium correlations and dissipative 
forces between two obstacles located far from each other, 
depending on the distance between them | r lk | and their 
mutual alignment 9i k with respect to the external field 
g. In view of (IB2QI) (IB22I) . the influence of the fcth ob¬ 
stacle on the fcth one is not equivalent to that of the fcth 
obstacle on the *th one (6 kl = 7 T — 9 ik ), i.e., these correla¬ 
tions are not reciprocal, 6n k i ^ Srii k , and the forces are 
non-Newtonian, f ki / -f ik . 

As is seen from (IB22I) , dissipative forces acting between 
inclusions are expressed, in the dipole approximation, in 
terms of induced density /r° of isolated inclusion. Distri¬ 
bution n°(r) = no + Sn°( r) for a single obstacle in 3D 
case takes the form 

„ r P q( r - r ')- 9 l r - r ' 

” (0 =n0 + J s 4 7r|r — r'| ^ (r>) ^ (B23) 

Far from the obstacle, when |r| a (a is its character¬ 
istic size), we can easily extract the leading asymptotics 
for the gas density perturbation Sn induced by the exter¬ 
nal held q = (1/2 — no)g: 

e qr-g|r| _ 

M r ) ~—O—7(r, q), (B24) 

[compare with Eq. mi Hr , a) is responsible for the 
sign of Sn dependence on r direction and, in turn, de¬ 
pends on |r| through a power law. 

Behavior of I(r,q) in the case of a spherical obsta¬ 
cle is defined by the asymptotics of the Bessel function 
K m+h {qr) (ZH: 

p~ qr(l —P cos 6) 

Sn° « V2na 2 -J(r, q), (B25) 


I ( r fcj, q, Xfc) = / 

J St 


(l+X fci Ufej)/ii(x , i ) 


+ 


e -?|x fe fc|+qx fcfc 

S k 4 7r|xfcfe | 


- A / W (Xfc , x •, r ki ) /i? (x')dx' fc 


dx'. 

(B21) 


As it follows from Eq. (IB21I) . I (r k i, q, x. k ) has a power- 
law dependence on r ki and in the case of a,; -C Ti k de¬ 
pends only on the mutual alignment of the obstacles with 
respect to the external Held g, i.e., on 6 ki , the angle be¬ 
tween r k i and g. 

Using expression (1B20I) for the density perturbation 
6n k , we can represent the force exerted on the fcth inclu¬ 
sion by the fcth one in the form that is similar to (IA13I) : 


fki 



n(x fc )<5?i fei (x fc )dx fc 


e -ij(l-/3cos0fci)|rfci| 


47 r|r fci | 



n(x fc )J (r fci ,q,x fc ) dx fc . 

(B22) 


E 


I ~ > atr, 

m —0 


1 + 


m 


m 


2qr 



7 m(COS 0 ), 


(B26) 

where a m = a m (qa) depends only on the obstacle 
radius a k and external Held q. The coefficients a m 
are from the Legendre polynomials expansion /r°(0) = 
e Pqacos6 a n P n (cos9) at the obstacle surface and 
can be obtained as a solution of equation (IB6I) . 9 is the 
angle between r and g and /? = (1/2 — no)/|l/2 — uq\ = 
±1. Distribution 5n°(r) for an isolated circular inclusion 
in 2D case was obtained in [45]. In the particular case 
of no < 1/2, the dipole approximation gives the follow¬ 
ing distributions for the gas perturbation: ahead of the 
obstacle, qr = — qr, 


Sn[r) w be 2qr j [3c + 1] (qr) = + [c + 1] (qr) *j, 

(B27) 

and behind it, qr = qr, 

Sn(r) « -fc|(qr) - ^ + ^ [2c + 1] {qr)~i |. (B28) 
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Figure 8. Asymptotic behavior of dissipative forces (a) /f 2 
(longitudinal alignment) and (b) /f 2 (transverse alignment) 
at large inter-obstacle separation ri2. The slope on (a) corre¬ 
sponds to the asymptotics /f 2 ~ r(" 2 3//2 . Equilibrium concen¬ 
tration no = 0.8, external field g (|g| = 0.5) is directed along 
the ir-axis, the impermeable circular obstacles are of radius 
a = 7 (in units of £), forces are in units of kT /£. 


Here, constants c = 2 [3(ga)A"i(ga)/ 2 (ga)] 1 and b = 
\/87mo(l - no)|l - 2no| -1 /2(ga)A' 2 " 1 (<7a) are expressed 
in terms of the modified Bessel functions I n and A' ra ; a 


is the radius of the impermeable obstacle (A —> 1). In 
the case of the point-like inclusion, qa ~ qi -C 1, this 
method gives Sn ~ e -29r r -1 / 2 for a region ahead of the 
inclusion and Sn ~ — r ~ 3 / 2 for the tail asymptotics. This 
is in qualitative agreement with the numerical results fill ] 
and coincides with the asymptotic behavior of the wake 
relaxation for a moving intruder [22| . The general form of 
the dissipative force in 2D case is analogous to Eq. (IB22I) : 

e -9|r fe i|+qr r 

fki ^-|7p2— J n (x-k)I(.r k i, q,x fc )dx fc . (B 29 ) 

It is easy to show that for longitudinal alignment 
ui 2 (ri 2 ) = —ri 2 /|ri 2 | 2 and the leading asymptotic be¬ 
havior fi 2 ~ that is in agreement with nu¬ 

merical result for nonlinear Eq. ©, see Fig. [SJa), when 
the distance between obstacles |ifoj| is much larger than 
their radii aq. The form-factor A depends only on exter¬ 
nal field g and the obstacle radius a. For the transverse 
alignment, the force leading asymptotics behaves expo¬ 
nentially, In ff 2 oc — <z|ri 2 | + • • •, that is also in qualitative 
agreement with numerical result, see Fig. [5Kb). 
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